Signatures of Strong Correlations in One-Dimensional Ultra-Cold Atomic Fermi Gases 
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Recent success in manipulating ultra-cold atomic systems allows to probe different strongly cor- 
related regimes in one-dimension. Regimes such as the (spin-coherent) Luttinger liquid and the 
spin-incoherent Luttinger liquid can be realized by tuning the inter-atomic interaction strength and 
trap parameters. We identify the noise correlations of density fluctuations as a robust observable 
(uniquely suitable in the context of trapped atomic gases) to discriminate between these two regimes. 
Finally, we address the prospects to realize and probe these phenomena experimentally using optical 
lattices. 



Strong correlations in low-dimcnsional Fermi systems 
give rise to interesting many-body states such as those 
responsible for High- Tc superconductivity and fractional 
quantum Hall effect (FQHE) in two dimensions (2D), and 
Luttinger liquid 0,(3 in one dimension (ID). In partic- 
ular, the Luttinger liquid phase has been the paradigm 
of low-energy physics in ID systems for about half a cen- 
tury ,3, 4] . This phase is characterized by the absence of 
fermionic quasi-particles even in the presence of a well de- 
fined Fermi surfaces with the relevant modes of dispersion 
represented by bosonic spin and charge excitations prop- 
agating at different velocities (spin-charge separation). 
Very recentlyja new regime, the spin-incoherent (SI) Lut- 
tinger liquid 0, Q , has become an active area of research. 
In contrast with the spin-coherent (SC) Luttinger liq- 
uid, here the spin-incoherence results from the induced 
spin-spin interactions scale J = 27r^?i'* (p)'^ / (3m^t7iD) 
being the smallest scale in the system, J <C T ^ Ep, 
where Ep is the Fermi energy, T the temperature, (p) 
the total density and Uid the interaction strength for 
contact interactions. In systems of ultra-cold atoms, this 
regime can be reached by increasing either the interac- 
tion strength among particles or by reducing the density. 
Thus, even at extremely low temperatures, fluctuations 
may drive the Luttinger liquid to lose spin coherence and 
only charge excitations remain as the dominant propa- 
gating mode. While there exists a growing experimental 
evidence in support for the SC Luttinger liquid [2], ex- 
perimental evidence for the SI Luttinger liquid is very 
scarce [1, Q . Even though it is possible to change den- 
sity in ID condensed matter systems by applying a gate 
voltage, lack of tunability in the interaction strength lim- 
its possible experimental realizations. In contrast, in 
trapped ultra-cold atomic systems the SI Luttinger liquid 
regime is unavoidable since low-densities are inevitably 
reached near the confinement edges. 

With the recent developments in the techniques of 
trapping and manipulation of ultra-cold atomic gases, the 
study of the strongly correlated regime of many-body sys- 
tems has acquired new momentum. Trapped atoms form 
an ideal many-body system that can be configured with 
extreme control and purity allowing for a thorough study 
of many-body properties previously inaccessible in solid 
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FIG. 1: Schematic sketch of the 3D cigar-shaped cloud (up- 
per) and array of ID tubes (lower) formed by optical trap 
lattice potentials. Yellow (light) regions represent the spin- 
coherent Luttinger liquid regime, while blue (dark) regions 
indicate its spin-incoherent counterpart. 



state configurations. Experimentalists are able to achieve 
strongly correlated regimes either by configuring the in- 
teraction to be strong by use of a tunable atom-atom 
scattering resonance, for example a magnetic Feshbach 
resonance fl^ , or by introducing degeneracy in the single 
particle ground state by creating a lattice trapping po- 
tential The former occurs when a two-body bound 
molecular states is made resonant with the open channel 
threshold by dialing-in an external magnetic field. This 
allows for the effective two-body interaction to have a 
negative (attractive) or positive (repulsive) sign with a 
strength that is dependent on the distance of the bound 
state from the scattering threshold. The lattice potential 
is induced by creating a standing wave pattern when two 
laser beams interfere at an angle [11]. The depth of the 
individual lattice well and the distance between wells is 
directly related to the intensity and the wavelength of the 
laser light. Such controllability made available by optical 
lattices has played an indispensable role in the trapped- 
dilute-gas demonstration of condensed matter phenom- 
ena such as superfiuid to Mott insulator transition 
In other experiments, multiple beams were used to cre- 
ate an optical lattice in two spatial dimensions, for ex- 
ample X and y, resulting in a configuration consisting of 
an x-y array of tubes as depicted in Fig. [TJ Each tube 
can act as an independent elongated quasi-lD trap if the 
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tunneling between the tubes is adjusted to be zero with 
the transverse trapping frequency large compared to all 
scales, including interaction energy scale. Such controlled 
ID-trap arrays allowed for the successful demonstration 
of the Tonks-Girardeau Gas regime l^, M\ and play a 
key role in theoretical proposals of achieving a dilute-gas 
SC Luttinger liquid 15j, |l6[. Moreover, trapped multi- 
species/spin atomic systems offer the added advantage of 
individual addressability. This allows for measurement of 
properties that are species/spin dependent. For example, 
it is possible to measure the spin up-up and down-down 
correlation function separately. 

In this Letter, we extend the ID Fermi gas studies to 
both regimes of the Luttinger liquid in two-component 
quasi- ID Fermi gas as confined in an array of tubes as 
shown in Fig. [T] Previous theoretical studies of this 
regime in ultra-cold Fermi gases have been limited to 
investigating the separation between charge and spin de- 
grees of freedom [Ta, ^li^]. Here, we focus on density fluc- 
tuations and show that there is a qualitative difference 
between the SC Luttinger liquid and its spin-incoherent 
counterpart. 

We consider a two-component ultra-cold atomic gas 
in a 3D magnetic trap. To define quasi-lD tubes, 
an optical lattice is imposed in the transverse direc- 
tion. The combined effect of the trap potential can 



be modeled by V{x) — 



Vo [sin^(27rx/A) -I- sin^(27r2//A)] . Here, uj±,lUz are the 
trap frequencies for the cigar-shaped {uJz < uj±) trap, A/2 
is the period of the optical lattice and Vo is its depth. 
The interaction between atoms is modeled by the con- 
tact pseudopotential Mnt(a?i —2^2) = U^dS^xi —X2), with 
C^3D — 4:nh?a3Y)/m, where is the 3D s-wave scatter- 
ing length. For large depth of the optical potential, the 
3D cloud splits up into ID tubes and each tube is then 
described by the following ID Hamiltonian: 



i'a(z) + 



UiD / dz%l)Uz)%lj\{z)i)i{z)ip-^{z), 



(1) 



where V'i(^) [V'trla?)] is the creation [annihilation] op- 
erator for particles of spin a. Vtiap(^) = mujlz'^/2 is 
the ID trap potential and C/id — —ifi^/maiu represents 
the effective ID interaction strength defined in terms of 
the effective ID scattering length oid = — (aj_/a3D){l — 
Ca3D/a±), with C = |C(1/2)|/V2 [13, [H- As we see, 
the ID interaction strength exhibits a "geometric" res- 
onance for a certain 3D scattering length, that can be 
employed to tune it by varying magnetic field. Here 
aj_ — {h/mujQ)^/'^ is the harmonic oscillator length in 
the transverse direction for the quasi- ID tubes with 
LjQ = (87r^Vb/7TiA^)^/^ . The number of particles in each 
tube is determined by the global trap imposed on the 3D 
cloud. To justify the above continuum description one 



should demand that all relevant energy scales of the sys- 
tem (collectively denoted by E), which can be tuned to 
a great extent, obey the inequality hjJz <^ E <^ hojQ. 

Now, we begin by discussing the homogeneous situa- 
tion i.e., Vtrap = 0. The model given by Eq. (P) (and 
its lattice variant, known as the Hubbard model) has 
been exactly solved by Bethe Ansatz [l^ and it is known 
that for repulsive interactions the low-energy and long- 
wavelength physics is described by the SC Luttinger liq- 
uid, with two copies of non-interacting charge and spin 
bosons, and is governed by the Hamiltonian [l| : 



Hi 



SC 



(2) 



which describes the fluctuations of the charge and spin 
densities above the ground state values. The bosonic 
field operators obey canonical commutation relations 

[ipaiz),Ila'{z')] = iSaa'S{z ~ z') and [ipa{z),ipa'{z')] = 

[H„(z),H„,(z')] = 0, where H„ = j^dt ipa is the mo- 
mentum conjugate to the bosonic field ipa and Wc(s) is 
charge (spin) velocity. For spin-rotation invariant sys- 
tems, the SC Luttinger liquid interaction parameter in 
the spin sector is i^s = 1 ' 1] . The spin and charge veloc- 
ities and the interaction parameter in the charge sector 
can be determined from the exact solution for any inter- 
action strength and density [20]. In the weak coupling 
regime, they can also be obtained by bosonization [ij. 
In contrast, the SI Luttinger liquid is described only by 
the charge bosonic field [2l|, \2^. The Hamiltonian in 
this regime, Hsi, can be easily obtained from Eq. ^ 
by dropping the spin part and making the substitutions: 
Kc = 1/2 and Vc — 2vp (which is exact for local interac- 
tions). We assumed that the temperature is very small 
compared to Fermi energy and considered the limit J — > 
first and then T — s- (opposite order of limits than in the 
SC Luttinger liquid). The behavior is equivalent to that 
of non-interactin g sp inless fermions in agreement with 
the exact results [lOj. 

We shall now include the effects of the trap on a quasi- 
lD tube via local density (Thomas-Fermi) approxima- 
tion. Thus, the average density of the system as a func- 
tion of coordinate can be found as a solution of the equa- 
tion 



dE[{p)] 
d{p) 



M - "^trap(2:), 



(3) 



where E[{p)] is the ground-state energy of the uniform 
system per unit length and fjL is the chemical potential 
fixed by normalization, J {p{z)) dz = N . In the SC Lut- 
tinger liquid regime, for weak interactions (near the cen- 
ter of the trap) [l^ 



Esc[{p)] 



h^n^ipy , UiDipY 



24m 



(4) 
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then from Eq. ([3]) the density reads 



where R = (2^/TOwf)i/2, (po) = {Sm^/h^T:^)^/^ is 
the density of the non-interacting uniform system and 
9 = Um/fiTrvFo, with vpo = hn (po) /2m. Due to 
the coordinate dependence of the density of particles, 
the spin and charge velocities, as well as the Luttinger 
liquid interaction parameters become coordinate depen- 
dent [15]. From bosonization results we can write that 
Vsiz) = vf{z) - UiD/2'Kh, Vciz) = vf{z)/Kc{z) and 
Kc{z) = {l + UiD/hTTVF{z))-^/'^. These resuhs are valid 
for weak coupling only, which is apparent since the ex- 
pressions break down near the trap edges where the in- 
teractions become stronger. Since the trap potential does 
not break spin rotation invariance, Ks{z) — 1. For strong 
interactions, the SC Luttinger liquid regime can only be 
realized near the trap center. Thus, we will approximate 
Kc{z) to its value at z = 0, KdO) = Kc.o; cf. Ref. [23 1. 

We can do similar calculations for the spin-incoherent 
regime with energy [l^ 



Esi[{p)] 



6m 



(6) 



since in this limit atoms behave as spinless particles. Us- 
ing Eq. ([3]) again, the density reads this time 



{p)si (z) = (po) \/l 



\z\<R, 



(7) 



where (po) = {2mfi/h'^n^y/^ is the density of the uni- 
form non-interacting spinless system. For the SI Lut- 
tinger liquid parameters we have: Kc{z) — 1/2 and 
Vc{z) = 2vf{z). Taking the above into account the 
Hamiltonians for the SC Luttinger liquid can be written 
as 



SC- 



Va,0 

2 



dZn 



1 



Ka,0 



(8) 



where dZa = dz/Va{z), Va{z) = Va{0)Va{z) = VaflVa{z), 
^Pa{Za) = ^Pa{z{Za)) and tla{Za) = VaiziZa))Ila{z{Za)) ■ 

One can easily check that the commutation relations of 
the new bose fields remain canonical. Hsi can be recov- 
ered from Eq. ^ proceeding as above. As we see, in 
the new coordinates the conformal symmetry is restored. 
We can therefore use this powerful property to calculate 
physical observables. Since the velocities vanish at the 
trap edges, no particle and spin current can flow out of 
the system and thus open boundary conditions (OBC) 
are effectively imposed on the system. 

Now we turn to calculate the matrix correlator of den- 
sity fluctuations which we shall show to differ qualita- 
tively for different regimes of the Luttinger liquid. The 



choice of this observable is motivated by the versatility 
of trapped atomic systems, which give us the unique op- 
portunity to independently address spin components and 
study correlators that are typicall y ch allenging to mea- 
sure in solid state configurations 2J, |25|. The matrix 



correlation function that we are interested in is given by 



Gcrcr'{z,z') = {S Pcr{z)S p^' {z')) , 



(9) 



where dpcr{z) = pa{z) - {pa{z)) = 5pc{z)/2 + a5ps{z) 
is the fluctuation of the density of a species. Here we 
consider only the smooth (non-oscillatory) part of the 
correlation function and rewrite the correlator in spin 
and charge (i.e., total particle) densities 

Gaa'{z, z') = i {5pc{z)5pc{z')) + (TCr' {5 Ps{z)5 p s{z')) , 

where we have used the property of spin-charge separa- 
tion to drop the correlations between the spin and charge 
densities. Using 5pc{z) = \/2/'Kdz'4>c{z) and 5ps{z) = 
^Jl/2TTdzfs{z), the correlators are straightforwardly cal- 
culated. This is achieved by using the conformal transfor- 
mation w' = (L/tt) hiw—iL/2 {L = 2za{R)), which maps 
the upper half complex plane {w ~ vr + iz; < z) (in the 
imaginary time formulation), with an OBC imposed at 
2: = 0, to a strip of width L {w' — vt' + iz'; —L/2 < z' < 
L/2) Introducing the functions A± = 

and R^ — 64:z'^{R)va{z)va{z'), the correlation functions 
for the SC Luttinger liquid are given by 



G^^,(z,z') = Gc{z,z') + aa'Gs(z,z'), 



where 



Ga{z, z') — ^ 



1 



1 



Rl Vsin^Aa cos^AJ 



(10) 



(11) 



For the SI Luttinger liquid the average in the spin sec- 
tor is zero for large relative coordinates (|z — z'| ^ ao, 
where ao is the average distance between atoms), since 
the spin excitations are exponentially damped 2l| and 
the correlator reads 



G^i,(z,z') = --^ 



1 

2^2 



1 



1 



sin^ Ac 



cos2 A+ 



(12) 



Let us discuss these functions in more detail. Remark- 
ably, the off-diagonal correlators (cr ^ a') show a qual- 
itative difference between the SC Luttinger liquid and 
its spin-incoherent counterpart; namely, G^^(z,z') > 0, 
while G^^(z,z') < 0. (Notice that for a noninteracting 
system G^i (z, z') - 0.) 

Having obtained the correlation function for a single 
tube we are in a position to describe the response from 
the array of tubes, which is more faithful to the experi- 
mental situation. Since we assume that there is no par- 
ticle transfer between the tubes, we can conclude that 
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FIG. 2: Qualitative behavior of G||/|Gc| as a function of 
the 3D scattering length, aao, in a harmonic trap for ^Li 
atoms {N = 10^,uj± = 21607rHz, loq = Utt x lO^Hz and 
Wz = 4707rHz) for the central tube (dashed) and the array of 
tubes (solid). As the interaction strength is increased, the re- 
sponse changes to negative, indicating the crossover from the 
spin-coherent Luttinger liquid regime to its spin-incoherent 
counterpart. Inset: "Phase diagram" of a ID ultra-cold 
atomic Fermi gas, the dashed curve is a crossover line be- 
tween the two regimes. 

they are uncorrelated. Thus, the many-tubes correlation 
function can simply be written as 

G%(z,z') = (Spl(z)5pUz')) = S,,G%{z,z'), (13) 

where the superscripts i and j label different tubes. 
The responses from different tubes simply add up. In 
the in situ experiment, when the measuring laser beam 
probes the whole array of tubes, one is able to measure 
the full integrated response, which is given by Gaa — 
J2i j!^iiG"^(z, z')dzdz'; this averaging provides an ob- 
servable with a good signal-to-noise ratio. In Fig. [2j we 
show the qualitative behavior of Gn/ \Gc\ for an array 
of tubes confined in a global harmonic trap. We have 
used realistic experimental parameters for ^Li atoms to 
see that it is possible to identify the response for the 
two regimes without pushing the values to extreme lim- 
its. We see that the SI Luttinger liquid regime is realized 
(as expected) for strong interactions and thus yields the 
negative response, while the positive signal for weak in- 
teractions is due to the SC Luttinger liquid correlations. 

In conclusion, we have proposed an experiment to mea- 
sure signatures of strong interactions in ID ultra-cold 
atomic systems. We have shown that correlations of den- 
sity fluctuations (an observable easily accessible in atomic 
physics experiments compared to their condensed matter 
counterparts) qualitatively distinguish between different 
strongly correlated regimes in these systems and there- 
fore provide an ideal probe for detecting these regimes 
in cold-atom experiments. While the off-diagonal cor- 
relators are positive for the spin-coherent Luttinger liq- 



uid, they are negative for its spin-incoherent counterpart. 
This result is robust, does not depend on the details of 
the trap and should be easy to identify in the experiment. 
This is also a novel proposal to measure properties of the 
spin-incoherent regime in a non-condensed matter sys- 
tem, which opens up new avenues for the study of Lut- 
tinger liquids. In particular, it would be interesting to 
extend these studies to the spin-imbalanced case. 
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